function y=max1(n,k) 
y=17.86*0.42*k/(0.8+0.42*k)*1.22*10^11/(1.22*10^11+n)*n*exp(-1.6)*(1-exp(-2/3*(0.8+0.42*k)))+22.99*k/(0.8+k)*1.22*10^11/(1.22*10^11+n)*n*exp(-0.28*k-2.4)/(1-exp(-2/3*k-0.8))*(1-exp(-2/3*(0.8+k)))
